Bouncing solutions in Rastall's theory with a barotropic fluid 
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Using Rastall's theory of gravity we reproduce some features of the effective Friedmann's 
equation emerging from loop quantum cosmology. We determine a class of bouncing cos- 
mological solutions and discuss about the possiblity of employing these models as effective 
descriptions of the full quantum theory. 
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I. INTRODUCTION 

The existence of an initial singularity in the standard cosmological model is one of the most 
important problems to be treated by a theory of gravity. This problem comes together with the 
existence of singularities in the spherical solutions of General Relativity (GR), connected with the 
notion of black holes. But, while in the black hole case the singularity is covered by an event 
horizon, being inaccessible to an external observer, the cosmological singularity does not have such 
kind of protection, and the emergence of the universe from a singular state renders the initial 
conditions of the universe quite arbitrary and, strictly speaking, unpredictable. The existence of 
singularities in General Relativity is related to the so-called energy conditions [1] that any ordinary 
matter should satisfy. 

Classically, the initial cosmological singularity can be avoided by introducing ad-hoc fluids that 
violate the energy conditions. This may even imply to consider a phase where the total energy 
density is negative (violation of the weak energy condition) or a negative pressure such that p+p < 
(violation of null energy condition). Another possibility is a modified theory of gravity, like scalar- 
tensor theories or non- linear Lagrangians. Quantum effects may lead to the appearance of negative 
energy states, as it happens in the Casimir effect jjj. Such fact gives the hope that the construction 
of a quantum theory of gravity may solve the singularities problem of General Relativity. 

rm 

One possible modification of the classical theory of gravity is given by Rastall's theory [3|, |4J|. In 
the latter the usual conservation law represented by the zero divergence of the energy-momentum 
tensor is given up; instead, the divergence of the energy-momentum tensor is related to the Ricci 
curvature scalar. The main argument in favour of such proposal is that the usual conservation laws 
connected with the divergence of the energy-momentum tensor are tested only in the Minkowski 
space-time (more specifically, in a weak field limit). Even if Rastall's theory is incomplete, in the 
sense that it lacks a Lagrangian formulation (at least in the context of the Riemaniann geometry), 
it has an interesting feature: it reproduces in a quite phenomenological way, but in an explicit 
covariant form, a distinguishing property of quantum effects in gravity systems, the violation 
of the classical conservation laws [5]. In this way, Rastall's theory may be viewed as a classical, 
phenomenological way to implement some concepts that are typical of a quantum theory of gravity. 
While Rastall's theory has already been studied in the context of the present evolution of the 



universe 



6l4l0l], its application to the primordial universe is still lacking. 
It is generally expected that quantum effects may lead to the avoidance of the initial singularity 
through, for example, a bounce phase which may connect the present expansion phase of the 
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universe to a primordial contracting phase. Such bouncing solutions appear in the context of 
string theory [111 . 1 1 21 ] and quantum cosmology 

Loop quantum cosmology is one of the most interesting proposal in the direction of applying 
quantum gravity to cosmology 



15l4l9l|. One of its main results is the modification of Friedmann's 



equation leading to, 



Pc 



where H is the Hubble parameter, p is the density, and p c is a constant characteristic density that 
defines a possible bounce. This last term is absent in the usual Friedmann's equation and formally 
it corresponds to the effect of negative energy densities. This structure must, in principle, lead 
to singularity-free model. A class of oscillating, non singular solutions has been determined, for 
example, in reference [201 ], 

In this paper we intend to study a Rastall's primordial universe model that can mimic the effects 
represented in Eq. ([I]). This can be achieved by considering a specific equation of state in the context 
of Rastall's theory. We show that bouncing solutions do emerge, but other possible scenarios 
appear. We will try to classify all possible solutions, by explicit integration of the equations, using 
a dynamical system approach. This work may be understood as a first step in verifying if Rastall's 
theory may lead to consistent cosmological scenarios, which can be tested against observations as it 
has been done with string effective action fl^ . More fundamentally, such analysis, complemented 
by a perturbative study, may shed light if indeed Rastall's theory can be considered as a classical, 
effective formulation of quantum cosmological models resulting from loop quantum cosmology, as 
we have freely speculated above. 

The paper is organized as follows. In Sec. HTJ we present Rastall's cosmology and introduce a 
suitable equation of state for the primordial phase. In Sec. IIIIl we find explicit solutions for the 
scale factor and the energy density, for some particular combinations of the parameters of theory 
which render the equations exactly solvable. In Sec. \IV\ we perform a dynamical system analysis 
in order to check the general behaviour of the solutions, even for the cases which are not exactly 
solvable. Finally, we present our conclusions in Sec. [V] 
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II. BASIC EQUATIONS 

RastalFs theory is based on the following coupled equations Q, 0] : 



7-1 T 



7-1 



(2) 
(3) 



i.e. Einstein equations properly modified in order to account for a non-conservation of the stress- 
energy tensor (recall that Bianchi identities are still valid). For 7 = 1, GR is recovered. 

Under the assumption of homogeneity and isotropy, the space-time is described by the 
Friedmann-Lemaitre-Robertson- Walker metric, 



ds 2 



dt 2 - a(tf 



dr z 



1 + kr 2 



+ r 2 (d6 2 + sin 2 



(4) 



where a(t) is the scale factor function which describes the evolution of the universe, and k = 0, ±1 
is the curvature parameter, indicating a flat, an open or a closed spatial section. Hereafter, we 
suppose a flat universe k = in order to simplify the calculations. Observation also supports a flat 



or nearly flat universe 
read 



21 



221 ] . Friedmann's and the conservation equation in Rastall's framework 



3— = 8ttG 



p + 3-(p + p) 
a 



7-1 



(p ~ 3p) 



(5) 
(6) 



where the dot denotes derivation with respect to the cosmic time t. 

In order to close the above system of equations, it is necessary to define an equation of state. 
Our main objective is to find bouncing solutions, a behaviour at the beginning of the history of 
the universe which is predicted by some quantum gravity theories, in particular by loop quantum 
cosmology. Therefore, hereafter we consider the following barotropic equation of state: 



P 



ap + uop . 



(7) 



If u) < 0, this equation of state allows to recover, in the context of Rastall's theory, therefore 
classically, the structure of the effective Friedmann's equation of loop quantum cosmology. We 
remark en passant that such equation of state has form similar to the so-called modified Chaplygin 



gas (MCG) models 



23|426|] 
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With Eq. (J7J), Rastall's equations ([21 [3]) assume the following form: 



a 2 



3^ = 8itG L - 1— ^ [(1 - 3a) p - 3ojp 2 ] \ , (8) 
p + 3- [(1 + a)p + up 2 ] = 7 — ^ [(1 - 3a)p - 6upp} , (9) 

Ob i— 



that may be rearranged to give 

;,2 



3- = 4vrG { [(3 - 7) + 3a( 7 - l)]p + 3u,( 7 - l)/> 2 } , (10) 
• = _> (1 + a + up) 

P o P [(3- 7 ) + 3a( 7 -l) + 6a;( 7 -l)p] ' 1 j 

If a ^ — 1, then equation (fTT|) may be integrated in order to give an implicit solution 

p A {l + a + up)% = Xa~ 6 , (12) 

where we have defined 

(3 - 7) + 3a( 7 - 1) B (77 - 9) + 3a( 7 - 1) 

1+q ' w 1+a ' 



and A is a positive integration constant. If a = —1, then Eq. (Ilip becomes 

- 3 ~ 27 + 3(7 - 1) In a = -31na + A, (14) 

LOp 

where A is another integration constant. Equations (112f) and (114j) can be analytically solved only 
for some specific cases, namely: 

(i) A=g -1,7=|; 

(«) A = -f 7 = 1; 

(m)a^-l, A = ^a^-1, 7 = fci); 

H a -1, § = ^ a ^ -1, 7 = ^ ; 
( w ) a = -1, 7= § ; 
(yi) a = —1, 7 = 1. 

In the following section we shall focus on the solution of each of the cases listed above. 

III. INTEGRATION OF THE EXACTLY SOLVABLE CASES 

A. Case 7 = 3/2 

In this instance we have A = B/uj = 3/2, and Eq. (|12|) becomes 

p(l + a + ujp) = Aa" 4 , (15) 
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whose solution is 

' = -^H 1+ (iT=H- <i6) 

With the above solution, Eq (fTOj) reduces to 

(ad) 2 = 2vrGA . (17) 

Taking the square root and integrating, one obtains 

aocti Hcxt" 1 , (18) 

which coincides with expansion behaviour of a radiation-dominated universe in standard cosmology. 
This result is quite curious, because such behaviour does not depend on a or u, which are free 
parameters. 

As for the sign in the parenthesis of Eq. (|16p . considering a > — 1, there are three possibilities: 
uj > with the minus sign; uj < with again the minus sign and u) < with the plus sign (u > 
with the plus sign would result in a negative energy density). In the first situation, the density 
decreases with time, tending to zero in the future. In the second situation, the density starts with 
a finite value and decreases to zero in the future. Finally, in the third situation, the density starts 
with a finite value and increases with time, tending to a value twice as big as the starting value. 
Considering a < — 1, then the density is non- negative only if uj > with the plus sign in the 
parenthesis of Eq. (16). In such situation, the density decreases with time, tending to a positive 
value in the future. 



B. u± -1, 7 = 1 

From Eq. (]13p one has that 



1 + a uj 1 + a 



and Eq. (|12p can be solved obtaining 

A(l + q) 



Since 7 = 1, Eq. ([10]) reduces to the standard Friedmann equation: 



(20) 



d 2 _ SvrGAjl + o) 
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If oj < 0, since the energy density must be positive, we must have a > — 1 which implies the 
existence of a singular state a — > 0, but with finite energy density p = {1 + a)/\oj\. On the 
other hand, if oj > 0, we have two disconnected evolutions: i) if a < —1, then a 3( - 1+a ) < wA. 
Therefore, we have an evolution from a singular a — > state with finite density p = |1 + a\/oj to 
a configuration with maximum scale factor a = (o;A) 1 / 3 ^ 1+a ^ and diverging density; ii) if a > — 1, 
then a 3 ( 1+a ) > ojX and we have an evolution from a minimum scale factor a = (wA) 1//3 ( 1+a ) and 
diverging density to an asymptotic Minkowski universe H — > for a — > oo. 



a j= -1, 7 



3(q-i) 

3q-1 



In this situation, we have that ^4 = and, from Eq. (I13p . -B/o; = 6(7 — 1). Substituting these 
values in Eq. (fT2|) . we get 

1 + a 



A 

p = —a 1- 1 



(22) 



and Friedmann equation ()10p becomes 



a 
a 



4vrG(7 - 1) 



CO 



Xa 7-1 - (1 + a ) 



(23) 



Observe that 7 — 1 and ui must have the same sign in order to avoid an imaginary Hubble parameter. 
In particular, 7 > 1 if a < 1/3. Remarkably, the above equation can be integrated to give 

7-1 



k exp 




4ttG(1 + a) 2 



t + 



A 



1 + a 



(24) 



w(7 - 1) 

where /c is a positive integration constant. It is important to realise that, since 7 = 3(a— 1) /(3a— 1) , 
here we cannot obtain the GR limit, since it would require a — > 00. There are two main cases: 

1. If q < 1/3, then 7 — 1 and oj are positive. In this case, the scale factor decreases with time 
and vice- versa the density increases with time. Note that if —1 < a < 1/3, in order to 
have a positive energy density we must have a < [A/(l + a)] 7-1 , which is impossible from 
Eq. ()24p . Therefore, the case —1 < a < 1/3 is excluded. On the other hand, if a < —1, the 
scale factor varies from infinity (t — > —00) to zero at a certain finite time. Correspondingly, 
the density is always positive and grows from p — >• |1 + a\/u) to infinity. This behaviour is 
somehow reverted with respect to the one of standard cosmology. 

2. If a > 1/3, then 7 — 1 and oj are negative. In this case, the scale factor increases with 
time and vice-versa the density decreases with time. Indeed, the former varies between zero 
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(t — > — oo) to a = [A/(l + a)] 7-1 [t — > +00), where the Hubble parameter vanishes and 
therefore represents an asymptotic Minkowski universe. Correspondingly, the density varies 
from p — > (1 + a)/|a/| to zero. 



D/ 1 3(a+3) 



Here we have B/oj = and, from Eq. (fl3j) . A = 6(7 — 1). Thus, Eq. ([12]) simplifies to 

p = Aa"^ 1 , (25) 
and Priedmann equation fjlQf) can be rewritten as 



(a^r _1 d) 2 = 4vrGA(7-l) 2(1 + a)a^ + a;A 



Hence, the scale factor may be expressed by: 



2vrGA(l + a) t 2 



L0\ 



7-1 



(26) 



(27) 



7-1 2(1 + a) _ 

Let us consider some possible scenarios. Again, the GR limit 7 — > 1 is not allowed here. 

1. If a < —7/3, then 7 — 1 < 0. If a; > 0, this corresponds to some sort of bounce linking 

two singularities, since the scale factor varies from zero (t — > —00) to zero (t — > +00) with a 

7-1 



maximum value a 



uiX 
211+qI 



for t = 0. Correspondingly, the energy density varies from 
zero (t — > —00) to zero (t — > +00), with a maximum value p = 2|1 + a|/w at t = 0. 

If a < —7/3, but w < 0, then the range —00 < t < +00 is no more allowed because it would 
correspond to an imaginary scale factor and density. There are two disconnected evolution: 
In the first t varies from minus infinity to the negative value which nullifies the expression 
between parenthesis in Eq. (|27p . the scale factor varies from zero to infinity and so the energy 
density. In the second t varies from the positive value which nullifies the expression between 
parenthesis in Eq. (|27p to infinity, the scale factor varies from infinity to zero and so the 
energy density. 

2. If —7/3 < a < —1, then 7 — 1 > 0, but 1 + a is negative. If u < there is no way out, the 
scale factor and the density become imaginary. Therefore, we can only consider u) > 0. The 
time t cannot go to ±00, but only varies between the roots (finite) of the expression between 



parenthesis in Eq. (|27|) . Therefore, the scale factor varies from zero to zero with a maximum 
for t = 0. Correspondingly, the energy density varies from infinity to infinity 



21+a 



with a minimum p = 2|1 + a\/co at t = 0. 
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3. If a > — 1, then 7 — 1 and 1 + a are positive. If u < 0, the scale factor varies from infinity 

7-1 

_2(l+a)_ 

density from zero to zero with a maximum value p = 2(1 + a) /\uj\. This is a proper bouncing 



(t — > —00) to infinity (t — > +00) with a minimum value a = 0(1+ ) an< ^ ^ e ener Sy 



cosmology, connecting two asymptotic Minkowski universes. If u > 0, we have again two 
disconnected evolutions. In the first t varies from minus infinity to the negative value which 
nullifies the expression between parenthesis in Eq. (j27[) . the scale factor varies from infinity 
to zero and the energy density from zero to infinity. In the second t varies from the positive 
value which nullifies the expression between parenthesis in Eq. (|27p to infinity, the scale 
factor varies from zero to infinity and the energy density from infinity to zero. 



Note that in order to have a bouncing solution it is necessary a very strict relation between the 
constants a (which comes from the equation of state) and 7 (which is a parameter of Rastalks 
theory). It can be interpreted as a dependence of the violation of the conservation law on the fluid 
properties. 

E. a = -1, 7 = § 

If 7 = 3/2, then the first term in Eq. (jTSjl disappears, and thus, it simplifies to 

p = \a- 2 . (28) 

It is interesting to notice that this equation may be seen as a particular case of Eq. (|25|) when 7 = | . 
In spite of that, there is no bouncing behaviour in the present situation, since the supposition that 
a/-l was crucial in order to obtain the result in the case (iv). This means that case (v) may be 
considered degenerate form of case (iv). 

It should be noticed that, from the results of cases (i) and (v), the scale factor behaves according 
to a oc 1 2 whenever 7 = 1, independently onaew (that means, independently on the behaviour 
of the fluid). This fact is even more curious because it coincides with the evolution of the scale 
factor in standard cosmology for a radiation-dominated universe (this could be explained by the 
fact that in both cases the trace of the stress-energy tensor is vanishing). However, the density will 
behave differently, depending on whether a/- 1 or a = - 1. 



10 



F. a = -1, 7 = 1 



In this case, Eq. (iMjl becomes 



1 



3uj In a — ojX 

and Eq. (|10p reduces to the case of General Relativity, 

8vrG 



d 2 

= 8irGp = 
a 3a; In a — ujX 



Solving this equation one obtains 



a = A exp 



t 2 



where here A = e A//3 . The expression for the density is 



3 (2irGuj 2 t 2 ) 



(29) 



(30) 



(31) 



(32) 



IV. DYNAMICAL SYSTEM ANALYSIS 



In the last section we considered some exact solutions of Rastall's cosmology, one of which 
leads to a bounce. In the present section we address a deeper analysis via a dynamical system 
investigation with the scope of checking the existence of other bouncing solutions for different 
choices of the parameters. 

Consider the following dynamical system made up of the acceleration equation plus the energy 
conservation equation in Rastall's cosmology: 

QH(p + p) 



P 



H 



(3 - 7 ) + 3( 7 - l)c 2 
7-1 



-4vrG 



P + 



(p ~ 3p) 



where c 2 = dp/ dp is the speed of sound. Note that Friedmann's equation 



3H 2 



8ttG 



P~ -^TT-(P- 3 P) 



(33) 
(34) 

(35) 



plays the role of a constraint for the dynamical system, indicating which solutions are physical. 

Since in our case we have a fluid of type p = ap + ujp 2 (which means that c 2 = a + 2up), then 
the dynamical system may be rewritten as 

6Hp(l + a + ujp) 



P 
H 



(3-7) + 3(7-l)(a + 2wp) ' 
7-1 



-AttGp 



a + cop + 



1 — 3a — 3up) 



-H 2 . 



(36) 
(37) 
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The critical points are characterized by the simultaneous conditions p = and H = 0. In our 
analysis we shall consider only critical points at finite distance: 



1. H = 0, 

2. H = , 



p = 0, 

1 2a + (7- 1)(1 -3a) 



P 



U 



37-5 



3. H 



8vrG(l + q)(3 - 2 7 ) 
3a; 



1 + a 



a; 



The second critical point obeys the constraint (|35[) if and only if it is equivalent to the first or to the 
third critical point. Therefore it is sufficient to consider the first and third critical points. Remark 
that, for a = — 1, the critical points coincide representing the Minkowski space-time. Linearising 
the dynamical system (I33"j) - (l34")) one obtains 

6p(l + a + up) 



Sp 



+ 



5H 



5H 

(3-7) + 3(7-l)(a + 2wp) 

36w(7- l)Hp(l + a + up) 
[(3-7) + 3(7- l)(a + 2up)] 2 ~ (3-7)+3( 7 -l)(a + 2a;p) 

7-1 



6H(1 + a + 2up 



>5p 



-3HSH - 4vrG 



a + 2w/j + 



'1 — 3a — 6up) 



Sp. 



(38) 
(39) 



At the critical point (p = 0, H = 0), the linearised dynamical system reduces to 



and at (p 



l+q 



if 



5p = 0, 
6H = -4irG 



87rG(l+q)(3-27) ■ 
3uj 



a + -t— — (1 — 3a) 



Sp 



(40) 
(41) 



Sp 



SH 



6(1 + a) 



(3- 7 )-3( 7 -l)(a + 2) 



8vrG(l + a)(3 - 2 7 ) 



3u 



Sp 



3i/ _ 8 , C (l +a) (3-2 7 )^ + 4rG 
3u 



7-1 

a + 2- (7 + 3a) 



Sp 



(42) 
(43) 



Looking for a bouncing solution, the critical points must form a diagram in the phase-space like the 
one in FIG. [H Moreover, both critical points must have H = 0. This is possible only if 7 = 3/2, 
consequently the linearised system reduces to 



Sp = 0, 

SH = -ttG(1 + a)5p 



(44) 
(45) 
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FIG. 1. Cartoon diagram of a bouncing solution. 



for the critical point (p = 0, H = 0) and 



5p = 



(46) 



5H = ttG(1 + a)8p 



(47) 



for the critical point (p = — ~nr->H = 0). It should be noticed that for 5p > 0, the point (p = 
0, H = 0) is an attractor and the point (p = —^-^,H = 0) is a repeller, and for Sp < the reverse 
occurs. For this reversed situation, the phase space diagram would look as the cartoon one in 
FIG. [2j The analysis of the dynamical system does not lead us to any new solutions with bounce. 



In this paper, we study a classical cosmological model based on Rastall's theory of gravity 
that can mimic quantum effects typical of loop quantum cosmology. Through the introduction 
of a convenient equation of state, see Eq. ([7]), we reproduce repulsive effects as the singularity is 
approached during the evolution. We also obtain some exact solutions, one of which exhibits a 
bouncing behaviour, described in the subsection D-3, of section III. In order to obtain such bouncing 
scenario, a connection between the equation of state parameter a and Rastall's parameter 7 must 
take place. Besides, some interesting radiation-like (even without radiation!) scenarios free of 





FIG. 2. Cartoon diagram of a reversed-bouncing solution. 



V. CONCLUSIONS 
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curvature singularities emerge. We finally perform a dynamical system analysis confirming these 
results. 

The bouncing solution here obtained can be studied perturbatively using the full covariant equa- 
tions. This may allow to establish some observational constraints on the model. More important, 
such perturbative analysis may help to verify if the similarity between Rastall's theory and loop 
quantum cosmology is just an accident or may contain a deeper meaning. This is an important 
aspect of the problem, since the basis of Rastall's theory contains some ingredients of quantum 
effects in curved space-time. 

It must be stressed the some perturbative studies in the context of loop quantum cosmology 



rave been already carried out 
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311 ] . But, most of them focused an inflationary behaviour and in 
29| the authors consider gravitational waves in a bouncing scenario. It will be crucial to confront 
the results of these references with a full perturbative analysis of the Rastall's models exhibited 
here. We hope to address this question in the future. 
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